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Abstract 



^ , We study the behaviour of a specific system of relativistic elasticity in its own gravitational 

' field: a static, spherically symmetric shell whose wall is of arbitrary thickness consisting of 

(•~*) ' hyperelastic material. We give the system of field equations and boundary conditions within 

the framework of the Einsteinian theory of gravity. Furthermore, we analize the situation in 
CO ' the Newtonian theory of gravity and obtain an existence result valid for small gravitational 

' constants and pointwise stability by using the implicit function theorem. If one replaces the 

elastic material with a fluid, one finds that stable states can not exist. 
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1 Introduction 

Let us take any static spherically symmetric shell consisting of an elastic material. Picture the 
KM \ shell to be a sphere from which a sphere of smaller radius is cut out around the center of mass of 

i the bigger sphere and replace the resulting hole by vacuum. The elastic material is assumed to be 

' isotropic and homogeneous. At first, ignore the own gravitational field of the elastic shell. When 

"switching on" the shell's self-gravitating field one will observe a deformation of the static elastic 
shell. In order to describe this process one can make use of either of the following two pictures: 
the spatial description (Euler picture) or the material description (Lagrange picture). In any of 
the pictures the undeformed body, B, is represented by a three-dimensional differentiable mani- 
fold that is endowed with a flat metric {body metric), Gab — B — 1,2, 3. The spacetime, M, 
a four-dimensional manifold, is equipped with the spacetime metric gap — a,[3 = Q, 1,2,3. We 
assume that there exists a natural state or relaxed state of the static elastic shell, i.e. one should 
understand this state as to be strain- and hence stressfree. In the spatial description we choose 
local coordinates y°' on M and let the motion of the body be given by a surjective mapping z — 
the so-called deformation map — from the spacetime onto the body so that = z^{y°') are 
coordinates on B. The field equations describing the behaviour of the static elastic shell viewed 
from the point of the spatial description are the Einstein field equations with elastic matter source 
on M. For the material description we have to assume that the spacetime has a foliation consist- 
ing of spacelike hypersurfaces. The motion of the body is given by the deformation map / from 
the undeformed body onto a submanifold in spacetime, namely /(B) x R. The tensor field gai3 
solves the Einstein vacuum field equations on M//(B) x M and the Einstein field equations with 
elastic matter source composed with the deformation map / on /(B) x M. For the static case, 
the deformation map of the spatial description is the inverse to the deformation map /" — 
a = 1, 2, 3 — of the material description. 

*E-mail address: christianelosertSgmx . at 
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From the experimental point of view it is clear that stable states for the above described config- 
uration exist for elastic materials but do not exist in the case of fluids. It would be of interest 
to derive this from the theory of elastomechanics. In [3] an existence proof for general static, 
self-gravitating elastic bodies in the Newtonian theory of gravity was given which includes the 
static elastic shell as a special case. Here we follow a different approach to derive an existence 
result and, in addition, explicitly write down the field equations and boundary conditions for the 
static elastic shell in both the Einsteinian and Newtonian theory. 

In the following section we will introduce some important quantities of elasticity that will simplify 
the tackling of our setting. We will make some crucial assumptions on the material under consid- 
eration. We obtain a nice expression for the energy-momentum tensor in terms of typical objects 
of elastomechanics and give the Einstein field equations and boundary conditions in the material 
description. In section 3 we consider our system of equations in the Newtonian limit. Still, the 
Newtonian field equation and boundary conditions cannot be solved explicitly for general elastic 
matter. Therefore, section 4 is concerned with an analytical approach to prove that for small 
gravitational constants and in the case of pointwise stability stable states of the static elastic shell 
exist. In the appendix, we solve the linearised system of Newtonian equations and illustrate the 
behaviour of the static elastic shell by means of three concrete examples. 



2 Field equations and boundary conditions in the Einsteinian 
theory of gravity 

2.1 Energy- momentum tensor 

We derive the field equations, namely the Einstein equations, from a Lagrangian principle in 
the spatial description. We aim to present the energy-momentum tensor using quantities within 
the framework of the theory of elastomechanics. Let us start with the general definition of the 
energy-momentum tensor: 

Definition. The energy-momentum tensor is defined by 

^ SC ^ 

where C is the Lagrangian density. 

In our case, the Lagrangian density is given by 

where n denotes the particle number density. The first term in brackets gives the energy density of 
the relaxed state and the second term, w, is called the stored energy function density. It is a very 
important quantity of elastomechanics and we will enlighten this in the few following paragraphs. 
In general, the stored energy function depends on the deformation map, its derivatives with 
respect to the coordinates on spacetime and the coordinates on spacetime: 

r A a. 

One can rewrite the energy-momentum tensor using the expression for the Lagrangian density 
such that one gets: 

tap = CUaUp -I- Sap, 

where is the J^-velocity — a timelike, future pointing vector field fulfilling u"§|^ = and 
gapu'^u^ — —1. The term Sap is called the stress tensor and satisfies Sapu^ = 0. The stored 
energy function determines the terminology of elastic materials, for example via the following 
definition: 
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Definition. // there exists a stored energy function such that the stress tensor takes the form 
then the material is said to be hyperelastic. 

All information about the specific classes of materials considered is contained in the stored energy 
function. Reminding of the form of the Lagrangian density one finds that we are actually dealing 
with hyperelastic material. Next, we want to introduce another important object of elasticity 
that is known as the Cauchy-Green strain tensor. Its name is justified when considering the 
fact that in spherical symmetry the {R,R) component of the strain tensor (H^b) measures the 
compression, {H^r) < 0, as well as the stretching, {H^r) > 0, of a material. In general, we 
have: 

Definition. The strain tensor is defined as 

■ ^ dy^ dy0 ' 

and 

Covariance of w under spatial diffeomorphism gives us 

Let us now assume that the stored energy function depends only on some invariants, J,;, of the 
strain tensor H^b'- 

w — w[Ji]. 

We make use of a certain choice of invariants^ of the strain tensor, H^b, namely: 

Ji=tr{H^B), 

J2^^[tr{H^Bf-tr(H^cH^B)], 
J3 = det(iJ'^i3) = 

At natural points (i/"\o,Y^\o, z^\o^ the strain tensor H"^^ takes the following expression 

in other words, the state of the static elastic shell at natural points is strainfree. In the relaxed 
state the chosen invariants of the strain tensor reduce to the following constants 

Ji\o = J2I0 — 3, J3I0 = 1- 

The stored energy function has to vanish and have a minimum in a locally relaxed state of the 
matter. From the expansion of the stored energy function at natural points we see that for 
isotropic, homogeneous materials there are constants A and /i — the Lame constants — such that 

dH^BQHCD = J^^^^^^C'D + 1^Hc(aHb)D. 

We reduce our analysis to isotropic and homogeneous material. Furthermore, we require our 
system to fulfill pointwise stability. In the case of isotropic, homogeneous elastic material — see, 
for example, |H] — this condition can be expressed as 

At > 0, 3A + 2/^ > 0. 



-'For the verification of the fact that the square root of the determinant of the strain tensor is equal to the 
inverse number density we refer to |4]. 
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Now, we will use all of the aforementioned assumptions and definitions to rewrite the energy- 
momentum tensor in terms of objects typical of the theory of elasticity. Due to the form of the 
boundary conditions it is more convenient to turn to the material description. Therefore, we 
derive the energy- momentum tensor in the material description, T"'^, from the one of the spatial 
description, tap^ and then compose it with the deformation map The energy-momentum 
tensor in the material description reads 



T"'3 = (J3 



dJi ' 'dJ2 dJs'dX'^dXB 



where the quantities Ji — i = 1,2, 3 — ,10, and Hab are to be understood as the analogous 
objects in the material description to the ones we treated earlier in this section in the spatial 
description. The symbol H^^ is meant to be the inverse of Hab- Note that H^^ is different to 
H"^^. The latter results of rising the indices of the strain tensor in the material description Hab 
with G^^. 



2.2 Field equations 

Having obtained a nice expression for the energy-momentum tensor in the material description we 
now turn our attention towards the equations describing the process of deformation of the static 
elastic shell taking into account the influence if the shell's self-gravitating field. The required field 
equations in the material description are the Einstein vacuum equations outside the deformed 
body and the Einstein field equations with an elastic matter source composed with the deformation 
map on the deformed body. That is, on M//(B) x M we have 

G"^(X^) = 0, 

and on /(B) x M 

where G"'^ is the Einstein curvature tensor, k :— Q is the gravitational constant, c is the 
speed of light and T"'^ is the energy-momentum tensor. 

We choose coordinates X"^ ~ (i?, 9, (f>) on the body B and introduce on spacetime M the co- 
ordinates a;" = (ci,a;") where = /"(X^) = (r = F{R),'d = 0,(p = (/)) — F{R) being a 
monotone function. The exterior Schwarzschild metric^ is matched to the outer boundary of the 
shell, that is where R = Ro is the outer radius. The hollow region in the centre of the shell (where 
R < Ri, Ri is the inner radius) is described by a flat metric. The body metric reads 

Gab = diag{l, R^, R^ sin^ 9). 

We divide the spacetime in three regions corresponding, respectively, to the hollow centre, the 
deformed body itself and the Schwarzschildean exterior of the shell: 

9 0.0 = diag{-C, 1, F^{R),F^{R) sin^ 0), 

9 af3^ diag{~A[F{R)], B[F{R)], F^{R), F^{R) sin^ 0), 

9 ^t3 = d^agi~ (^1 - , (^1 - , F^R), F\R) sin^ 0), 

where G is a positive constant and M is the central mass. 

We require asymptotic flatness, that is 

gtt[F{R)]^-l asF^oo, 
grr[F{R)] — * 1 as F ^ 00. 

■^This makes sense since we will require asymptotic flatness and we, therefore, take into account the Birkhoff's 
theorem — see 
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In a slight abuse of notion we write T^'' =ffa^ T^^p, and get 



rp tt 



'dF y/BF"^ 
dR i?2 



9w 



f2 aw 



,dF 



F^ dw 
F^' 



dw 



i?2 / a J, 



,dF.2 F^ dw 



Keeping these expressions in mind, and after having computed the Einstein curvature tensor we 
write the Einstein equations on /(B) x R in the following way: 



,,T * — I (U 



dA 



1 1 



'^^''^ ~ ABF F^^BF^' 

dB OA &\A (dA_\'^ dA dB 

~ ^ ~ 2B^F^ 2ABF^ 2AB AA^B AAB^ ' 

(3) 

If the first and second Einstein equations — and — hold then one can, instead of equation 
0, consider — see, for example, — the only part of the conservation law that is not fulfilled 
identically, namely: 

V^Tr" = 0, 

which can be written as 

r fF(l)'^ , .dB dF 2 d^F ff^ 45^ AF dw dF d'w 

L A dF^dR' dR^ 2B R R^^dJi dRdJidR 

L i?2^ A ^ dF^dR> dR^ 2B i?2 R^hj2 



i?2 dJ2dR 

^ ,0BdF2 d^ |f|| Ajj^fBF^ _ 4B^ dw 
Li?4^ A ^ dF^dR> ^ dR^ 2B i?4 R5 J^Ja 

Tails -'"-'^-'S-"^ w 

For convenience, we choose equation (0} to substitute for equation (|3Jl in our further investigations 
and consider equation together with the first and second Einstein field equation, and ||2J), 
the demanded system of field equations. 



2.3 Boundary conditions 

Clearly, additional equations will have to be fulfilled at the inner and outer boundaries of the 

deformed body in order to satisfy the standard matching conditions — see, for example, [7]. The 

fundamental boundary conditions of our system have to be derived from the condition that the 

1 

first fundamental form, ga/s, and the second fundamental form, Kafs, of the three metrics dat, 

2 3 

9ab, 9ab {o,,b — 1,2,3 Spatial indices) have to coincide at the inner and outer boundaries of the 
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static elastic shell. We consider the hypersurface TC : F{R) = const and the normal vector field 
of n", with 

ax" ' dr- 

The induced metric hap on Ti. is Lorentzianlike, Ha being spatial. 
Therefore, we have 

Kpdx°'dxl^ = gttc^dt^ + F{Rf{de'^ + sin^ edct)"^), 

where ^ denotes the Lie derivative. Considering the hypersurface Hi : F{R) — F{Ri) it follows 
that 

A[F{Ri)] = constant, (5) 
B[F{R,)] = 1. (6) 

The boundary conditions derived from the matching conditions at the hypersurface Tio ■ F(R) — 
F{Ro) read 

S|F(K.)1 = AIF{R,)]-K (8) 

As a consequence of the 3+1 decomposition of the manifold we also have to look on the constraint 
equations: 

^""^Rlh] + {K\f - KabK'^" - 2Tabn''n\ (9) 
D\Kab - K^hab) = {(j^n^'aiTdbn"), (10) 

where is the Ricci scalar with respect to the three-dimensional metric h, D'^Mab '■— 

h'^d'ha"' hb^ V** Ma'b' is the projective derivative, and {(t)*)°'b is the pull-back of the embedding. 

The momentum constraint equations (|10(l are fulfilled identically, whereas the Hamilton con- 
straint equation ^ leads to the following boundary condition: 

which is equivalent to 

^/|^r„n«|a„ = 0, (11) 

where n^ is the normal vector of the hypersurface R = constant. 

Henceforth, we are able to give the entire system of field equations — iQ), Q and Q — and 
boundary conditions — ©, O and — of the static elastic shell in its own gravitational field 
within the framework of the Einsteinian theory of gravity. Let us now investigate the situation 
in the Newtonian theory of gravity. 



3 Newtonian limit 

Our aim is to derive the system of the Newtonian equations from the results we obtained in 
the latter section. From the first Einstein equation ^ we derive the following form of the field 
B[F[R)\, namely: 

5[^(i^)].(l-^^™[^(^)]^" 

where 



c2F(i?) ; ' 

A-K r" - 

m[F{R)] = — / F\R){p^c^ + w[F{R)])dR. 

C JR, 
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We insert the latter expression for the unknown field B[F{R)] as well as the following form of the 
unknown field A[F{R)], namely 



,_2U[F{R)] 



A[F{R)] = c-^ exp 

— where U can be viewed as a potential — into the system of field equations (0), (01) 
and boundary conditions In order to obtain the system of field equations and boundary 

conditions in the Newtonian theory of gravity we take the limit c — *■ oo. In this Newtonian limit 
the first two Einstein equations ^ and Q reduce to the Poisson equation 

AU ^Ang{p)nerot, (12) 

and the conservation law Q in the Newtonian limit gives the load equation 

^ ((^^ ^newt {F^ )newt) ~t" ( ^~dR^^ ^newt — q-^ {p^newtj (-^^) 



where 



{Tr )newt — '^F^^dR^ 



dw dw F'^ dw 



(F^ ^newt — iF:p^\ncwt — ( j n ) 

art 

{P)newt = -F2{-T^)P0- 



dw ((dF.2 , F^\ dw , F^ ,dF^2dw 



dJi y-dR' R^)dJ2 R^^dR' dJj. 



F'^^dR' 

The remaining boundary conditions in the Newtonian limit are 

J^{Fr^)newt \dM= 0. (14) 

Equations 1)121) . (|13|l and H14|l form the complete system of Newtonian equations for the static 
elastic shell in its self-gravitating field. 

4 Main Theorem 

The resulting equations of the latter section — p2|) . (|13|l . p4|l — are too complicated to be 
solved explicitly for general elastic materials. Therefore, we will resort to analytical methods in 
the sequel. In order to use the machinery of the implicit function theorem we write our field 
equation and boundary conditions as a map between Sobolev spaces. Integrating the Poisson 
equation (|12() and inserting it into the load equation 113|) we find the resulting Newtonian field 
equation to be an integro-differential equation of the following form: 

■^{{F/)newt ~ {T-g^)newt) + (^) ^ {^^F/ ) newt = ^4 (^) ^ F^ {R)dR. 

The latter together with the boundary condition 

-^{F/)newt |aB= 

leads to modelling the characteristic mapping, T , corresponding to the static elastic shell in the 
context of the Newtonian material description. It is obtained from the system of field equation 
and boundary conditions in the Newtonian theory of gravity, that is 

Definition. S(tatic)E(lastic)S(hell) map: 

J" : M/2'2((i?,,i?„) X M) ^ W^'^{{R,,Ro) x R) x W^'^{{Ri,Ro} x M) 
[F{R), g] ^ T[F{Rl g] = {E[F{R)] - ge[F(i?)], 6[F(i?)]) , 
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where 



Static Elasticity Operator 

~. , . d^F 4g AF.dw dF d^w Ate^F^ 4(g) V 

E[Fm ■- [2— + _-_]_+ + [— ^ + — ^ - 

i?4 J dJ2 i?2 £)J25iJ ^ L ^4 ^4 

Force Operator 

e[F{R)] -.^^ j%\R)dR, 
Boundary Operator 

fr^.„M ^dF,dw ^F"^ dw F^ dw , , 
■■='dR^dT,-''R^8j:2'-R-^8Tj 

In the above definition W^''' stands for the standard Sobolev spaces of functions having p weak 
derivatives in L'' — see, for example, ^2^. 

Remark. The Sobolev spaces were chosen out of technical reasons. The first factor in the range 
of J- rises from considerations concerning the field equation, the second factor in the range from 
considerations concerning the boundary conditions. The weight associated with the weak deriva- 
tive in the first factor of the range comes from the appearance of second derivatives in the static 
elasticity operator while that of the second factor is dictated by a lemma in lOj. We will make 
use of this lemma later. 

Before giving our main result we will recall one of our crucial assumptions and we will consider 
some technical lemmata that will be used to prove our main theorem. 

Assumption 1. We require our system to satisfy pointwise stability. For isotropic linear elas- 
ticity pointwise stability — see, for example — is fulfilled when 

^>0 3/v = 3A + 2/x > 0, 

where A and fi denote the Lame constants. 

Remark. The constant k is called the bulk modulus. The inequality 3k > allows negative A 
with — |/i as a lower bound (auxetic materials) . 

Lemma 1. T e (^U,W°''^{{Ri, Ro) x M)) x W^•'^{{R^,Ro} x M)^, where U is an open subset 
of the Sobolev space W'^'^{{R^, Ro) x M). 

Proof. To check if the SES map is , we consider the special case of m = and p = 2 in the 
assumptions made in a lemma in ^HI- We find that one easily derives the validity of the assertion 
for the pair of the following operators, namely, the static elasticity operator and the boundary 
operator. To show that the force operator is C^, we make use of a corollary in ^ that reads as 
follows 

Corollary 1. If T : U C W^'^{{Ri,Ro) x R) -> W°'^{{R^,Ro) x M) x W^'^{{R^,Ro} x R) is 
-Gateaux then it is and the two derivatives coincide. 

Essentially, the force operator corresponds to -^U[F{R)]. The computation of the Gateaux 
derivative of ^U[F{R)] gives us 

l^{±,U[F^r){R)])r=, = D{^U[F^,){R)]) ■ X, 
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where 

= i? + rx(i?); Fo{R)=R 

The Gateaux derivative of ■^U[F{R)] reads 

l{^U[F,^,iR)])^^o = ^(1 + 2§)x = Di^U[F,om) ' X- 

It is an element of the space of hnear continous maps C{W'^''^{{Ri, Ro) x M), W°''^{{Ri, Ro) x M)). 
Now, R I— > Z)(^C/[F(o)]) is continuous for R e [i?i,-Ro] since the linear operator defined by the 
directional derivative is clearly bounded. 

Thus, the force operator is C^-Gateaux, from where it follows that the force operator e[F{R)] is 
C^. Therefore, we derive the SES map J^[F{R),g] — as defined above — is C^. □ 

Lemma 2. J^[Fo = R,go = 0] = (0,0). 

Proof. For J^[F = Fo,g = Go] = (0, 0) wo find Fo aiid Gn to be such that Fa = R and = 
as — inter alia — the invariants of {H^B)newt on which w depends are constant for Fo = R and 
Go = 0. We have 

J^[Fo = R,Go = 0] = 0- □ 

Lemma 3. DfJ^[Fo,Go] is an isomorphism from W'^'^{{Ri, Ro) x K) onto VF^'^ ({Ri, Ro) x M)) x 
VF5.2({i?.,i?„} xM). 

Proof. We divide the proof of this lemma in three parts: 

First task. From the Taylor expansion of — which is supposed to be at least C"""*"^ — around 
R for 6F sufficiently small we see that 

T[Fo + SF, Go] = DfHFo = R,Go = 0)SF + o(| 

The components of the linearised SES map read 



E [6F] = 



e [SF] 



b [6F] 



- fj^rs:T:^\ , ^ d ^^^^ 2{SF) 



(A + 2m), 



(i?^ - Rf), 



±{SF){X + 2t,) + ^X{6F) 



lin , 



lin 



hence JT [5F,G] := [E [5F] -G e [SF], b [SF]). 

Neglecting higher orders, wc are now able to concentrate on J^[Fo + SF, Go] to see if DJ^(Fq 

Go = 0) is an isomorphism. 
Second task. Next, we will show that T[Fq + SF, Go] is injective. 

lin 

We compute the solution of E [SF] = 0. The latter reads 

which is an ordinary differential equation of the Eulerian type. 
The solution of the above equation reads, 



SF{R) = uiR + U2R~ 
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where Ufc G M for k — 1,2. Inserting the solution SF in the hnearised boundary conditions 
we derive 

ui(3A + 2^) - ""2-^1 |9B = 0. 
We can write the latter equations in matricial form as 

Au = 0, (15) 

where 

The determinant of the coefficients matrix of this system of linear equations, det A, is 

detA= 12/t^(-i3 - -^), 

where 3k = 3A + 2/i. 

Assume for the moment that the above system (|15|l has solutions other than the trivial solution, 
that is equivalent to the fact that det A = 0. If this is the case it follows from pointwise stability 
that Ri — Ro- However, Ri = Ro is not an allowed solution. Thus, det A ^ and the system 

lin 

(|15|l has only the trivial solution. Therefore, [5F\ = has no kernel for Q ~ except for the 

lin 

trivial solution. That is, the mapping [SF] is injective for Q = Q. 

lin 

Third task. Now, we will show that [SF] is surjective for Q = 0. 
Let us write 

and suppose that 

e X wi^^{{R,,Ro}). 

The solution {SF)g{R) of the above differential equation, H16|l . consists of the sum over the 
solution {5F)h{R) of the homogenous differential equation, that is 4>{R) — 0, and a particular 
solution {5F)p{R) of the inhomogenous equation. The solution to the homogenous case as already 
computed before is 

iSF)h{R) ^uiR + U2R-^. 
One way of obtaining a particular solution to H16|l is by means of the following ansatz 

iSF)p{R) = v{R)R. 

After a few calculations we get 

Since we are searching for any particular solution, it is possible to choose the constants of inte- 
gration in a way that the last term in the above equation vanishes. Therefore, the solution of 
equation fT?)| reads 

{6F)g{R)=uiR + u2R-^ + Rj^(^J R <f)(R) dR^dR. (18) 
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Inserting this solution into the boundary conditions 



^i5F),{R) + l{5F),{R)]\aM^0, 



where 7 = 2A(A + 2/i) ^, we can determine the form of the constants ui and U2 in terms of v{Ri) 
and v{Ro): 



U2 



R^R^ 



(7-2)(i?3_i?3) 



ui = - (^«2(7 - 2)i?r' + ^^^"(^0 + 7«(^'^) ) • (19) 



Let us write 
where 



{SF),{R) = {SF),,{R) + {SF),.AR) + {SF),,{R), 

{SF)g,{R) = ui{(biR^),q}{Ro))R, 
{5F)g,{R)^u2{^{Ri),4>{Ro))R"^. 
{SF)g,{R) = v{cj,{R))R. 

Using the fundamental theorem of calculus it is easy to see that 

idF)g,{R) eW^'^{{R,,Ro)), and iSF)g,iR) e W^^^ {{R,, Ro)) . 

Concerning {SF)g^{R), the following is true 

A.iSF)gAR)(^C"{{R.,Ro)), 



J^"\^{SF)g,{R)\^dRy 



The same arguments for the second weak derivative of {SF)g^ (R) lead to the result that 

{SF)g,{R)£W^-'{{R.„Ro)). 

We conclude that the solution {5F)g{R) of the inhomogenous differential equation with 
4>{R) & W°''^(^{Ri,Ro)) X W^''^(^{Ri,Ro}) is an element of the Sobolev space as required, namely 

{SF)g{R) e W^'^{{R,,Ro)). 

lin 

Hence, [SF] is surjective for Q = 0. 

Remark. At this point let us say something about the possibility of taking the limit i?,; — > 0, 
which yields a static elastic sphere. It can be checked that 

= OiR-P), 

where p < 1, € W°'^{{Ri, Ro)) x W^'^{{Ri,Ro}). Therefore we conclude from ^ that 

= 3-p = 



v{R) = J R'^(^J R ' dR 

1 /R^-P R^P R}-p R}-p 



A~p\l-p 3i?3 l-p 
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and 

v{Ri)^0{Rl), RARi) = 0{RD, 
where g > 0. Thus, it follows from (|19|l that 

U2 — > as Ri 0, 

and the constant ui is bounded for Ri — > 0. The solution of the linearised system of equations 
for the static elastic sphere reads 

{SF){R)'^'="ui{R,^0)R + rJ^ ^"*(/ 'l>(R)^dR. 

lin 

Again, the linearised mapping [^F] for ^ = is both injective and surjective. Therefore, 

DFJ'iFa,ga) : W^-^{R,, Ro),R) ^ W°^^ {{R^, Ro),R)) x W^'^{{R^,Ro},M.) 
is an isomorphism. □ 

Now we are ready to state and prove our main result. 
Theorem. Let 

T : W'^'^{{R,,Ro) X K) ^ W°''^ {{Ri, Ro) x M) x W^'^{{Ri,Ro} x M) 
[FiR), g] ^ T[F{R), g] = {E[F{R)] ~ c;e[F(i?)], , 

where 

Eim)] - + _-_]_+ + [— ^ + — ^ - 

fr^/„M ^dF,dw ^F"^ dw F^ dw , , 

T/ien f/iere exists a neighbourhood Ng of g^, Ng C M and a neighbourhood Np of Fq, Np C 
W'^''^{{Ri,Ro) X R) and a map F G C^{Ng,NF) such that 

(i) T[F{R,g),g] = {Q,Q) ygeNg, 

(11) T[F{R), g] = (0, 0), [F{R),g] eNgX Nf, implies F{R) = F{R, g), 
(ill) F'{R,g) ^ -{DFT[F{R,g),g]y^ o DgT[F{R,g),g], wheregeNg. 

Proof. The proof of the above theorem is a direct consequence from Lemma 1, Lemma 2 and 
Lemma 3 as well as from the implicit function theorem — see n 
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5 Concluding remarks 



From the above theorem it is now clear that as long as pointwise stability is fulfilled, for example 
perfect fluids will not satisfy this condition, stable states of the static elastic shell in its own 
gravitational field exist if the body is sufficiently small. 

The present work can be extended in several directions. It would be of great interest to provide 
the analogous analytic proof for the existence of stable states in the theory of Einsteinian gravity. 
Furthermore, one could undertake some numerical investigations picking out some specific realistic 
material. It would also be possible to replace the vacuum inside the hollow centre of the shell by 
some matter, for example air. 
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A Solution to the system of linearised equations 

The Newtonian system of field equation and boundary conditions is too complicated to be solved 
for general material, still the linearised equations can easily be solved analytically. In addition, 
we will list a few examples for which we investigate the behaviour of the static elastic shell or 
the static elastic sphere as a special case taking into consideration various models of matter. We 
consider the system of linearised equations 



^iSF) + liSF)) =0, 
where 

IC = g-p^, and 7- 



R'^" 3(A + 2/x)' ' X + 2n' 

As we know from section 4, the solution to the above system for a general — see equations 
(dl), 113 and lO— reads 



{SF){R) = i^ui{cj,{R,),<j)iRo)) +v{cj)iR)) ] R + U2{^{R,) , (t>{Ro)) R^"" , (20) 

where 

« 1 / /■« =3 = = \ 



v{R)^ I ^{ I R m) dR ]dR, (21) 



U2 



R^R^ 



(7-2)(i?3-i?3) 



{R^^'"iR^) - ^o^v{Ro) + (1 + j){v{R^) - «(i?o))), 



(22) 

ui = - (^"2(7 - 2)i?r' + R^^^iR^) + 7«(^^)) • (23) 

Inserting (/)(i?) =IC{R- ^) into (EB we obtain 

R^ Rf Rf i?f 



10 2R 10i?3 2 
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And the following holds, 



\ irf^o ^ 

^«(i?o) = /C(^--^ + ^j. 



With these results it is clear that 



U2 = 1^-, „x , „n fyTT „ -Kq -t; — fTT T — -K, 



(7-2)(i?3-i?3) V 10 ° 10 i?3 2 ' 2i?„ 



(l+7)(i?3-i?f) V 10 ° 10 i?3 2 ' 2i?o 

After some simplifying computations we have, in the end, the solution to the system of linearised 
equations, namely 



where 



Eg . 5A + 6At 2 3(5A + 2/x) i?f Rf, 
'3K(i?3-i?3)i 10 ^° 10 Rl^ Ro' 



_ RlRl , 5A + 6m ^2 3(5A + 2a*) 2 i?3 

= 4Mi?g-i??) To — + ^ i?:^ ' 

and 3k = 3A + 2^. Note that for Ri — Ro, ui — U2 — 0- 

Next, we want to give some common examples of elastic materials and investigate in which 
regions the material is stretched and where it is compressed. The (R,R) component of the strain 
tensor (H^b) measures the compression, (H^jf) < 0, as well as the stretching, [H^ji) > 0, of 
the material. For the linearised system of equations of the static elastic shell, we have 

(H n)iR) = ^iSF)iR) = ^^^i—+U.-2U2R-^). (24) 

Example 1: Ideal cork. We consider a material where the Lame constant A = 0. This is a 
model for ideal cork'^. For such a material 

^i(A = 0) = -Yo^l\^ < 0, ^2(A ^ 0) = > 0, 

where Ri — bRo, b G [0, 1). If 6 = we have the special case of a static elastic sphere — see 
example 3. We apply Descartes' rule of signs'* to the following polynomial: 

,10 o 20 

and see that in this case we have at most one positive zero. Next, we consider {H R){Ri) and 
{H r){Ro) and investigate their signs: 

.,^.x=o2 gnplRl b'-l .=,2 g^plRlb\9-l) 



''See, for example, the following homepage: http : //home . att . net/ ~ numericana/ answer/physics .htm . 

■'Descartes' rule of signs states that for a given polynomial the number of sign changes of the cocffients of the 
polynomial gives the maximum number of positive roots — see e.g. [2). 
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Thus, we know that for this particular material, there exists no positive zeros and the material 
is compressed in the whole shell. Alternatively to the latter discussion, we can investigate what 
happens to the solution (SF){R) at the inner radius, Ri, and at the outer radius Ro- 

(6Fm) + ^i(A - 0)i?, + fi2(A 0)i?-2) 



10(1-63) 



- -i-Wb' + nb' + 66^ - 36^ - 4) < 0. 



10(1-53)' 

That means that under the influence of the gravitational field of the shell itself, the inner radius 
increases and the outer radius decreases. Therefore, the body is compressed. 

Example 2: Ideal rubber. Here, we want to investigate the model for ideal rubber^, where 
^ ^ 1. We have 

u,{fi « 1) = --^J^i-Sb" + 2b^ + 1) < 0, 
XR^ 

U2{fi < 1) = ~ 8^(1 -63) ^^^^^^ - + 1) < if ^ ^ 

U2{ti < 1) - -7-Tr^W&^(263 - 36^ + 1) > if A > 0. 
8^(1 — 6-^) 

From 

-|;^(2&^ - 36^ + 1) < 0, 

(/.)(i^o)^^<^-g^(263 - 36^ + 1X0, 
we see that for A > the material is compressed in the whole shell. 

Example 3: Elastic sphere. Specialising on the elastic sphere^, we let Ri ^ and see 

that 

?Ii(i?, ^0) = --^(l + -^) <0, and U2{R^^0)=0, 

where c^ — X~\- 2/i. 

Remark. The constant C2 is the speed of the propagation of the progressive, strongly elliptic 
wave — see, for example, 8 . 

Example 3 has exactly one positive zero. In the case of the sphere the solution reduces to 

47rrt^ / , T , 2 Co 



(^i^)«,^o(i?)^^^U3-i?i?^(l, 3^ 



^2 

Considering the latter equation, we see that within a sphere of radius R — Ro + |^) the 
material is compressed since 

{H r) < 0, 
whereas outside this sphere the material is stretched: 

{H r) > 0. 

^See http : //home . at t .net/j ~ numericana/ answer/physics. htm. 

^This special case has already been considered in and l^- 
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Note that if A = (ideal cork) we have the following 

{H r){Ro)Vo. 
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